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Abstract. Based on a unique waveform with strong exponential localization 
(— I ' property, an exact mathematical method for solving problems in signal analysis in 

O I- time-frequency domain is presented. An analogue of the Gabor frame exposes the 

^ fin non-commutative geometry of the time-frequency plane. Signals are visualized using 

■ graphical representation constructed. 
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1. Introduction 

Time-frequency (TF) representation of signals is essential in modern theories. A review 
■ of TF analysis methods can be found in pQ, for a brief introduction see [2]. Among special 
Q\ ■ mathematical methods of TF analysis a central role is played by Gabor analysis (see, 
e.g. [3]). However, there are difficulties regarding the practical implementations: (i) the 
existing methods of TF analysis are computationally complex, and (ii) understanding 
of the TF plane geometry is still poor. It is the aim of present paper to fill this gap. 

In TF approach to the sampling, the effects of localization are critical, especially, 
close to the quantum uncertainty limit. Localization properties of orthonormal bases, 
when viewed in global context of the whole of the TF plane, are subject to the restrictions 
imposed by Balian-Low theorem [I]. The theorem, as another form of uncertainty 
relations, reveals deep features of Fourier transformation, in its connection to the 
noncommutativity of the TF plane. 

The fact that this phenomenon has a very general character is known for decades. 
Parallels with phase-space quantum mechanics and quantum Hall effect describe the 
phenomenon best (see, e.g. [5]). As it is always the case with similar problems in 
different research areas, they pose a common mathematical problem too. Perhaps the 
best example in this course is the fcg-representation constructed by Zak [6], originally 
for the purposes of the solid-state physics; it was then developed into a general 
mathematical method used widely in the theory of signals. The first mathematical 
analysis along these lines goes back to von Neumann [7j, within the context of the 
phase-space quantum mechanics; the idea of elementary signals was put forward by 
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Gabor [5]. Balian [5] proved that the bases associated with the Gabor-von Neumann 
lattice are singular when the orthogonality additional condition is imposed. The similar 
fact in connection with the quantum Hall effect was established by Low [TO] . A way 
of circumventing the difficulties in realization of von Neumann's and Gabor's ideas 
posed by the Balian-Low theorem on localization was shown in [TT], by constructing 
an irreducible representation of the anticommutative subgroup of the Weyl-Heisenberg 
group, the group of translations-modulations. A unique entire function that naturally 
arises in that approach leads to exponentially localized basis which solves the problem 
of local asymptotic approximations to entire functions, in Hilbert space metrics [T2] . 
Due to its unitary invariant form, the approach can be applied directly to the signal 
analysis in TF plane. 

Conventionally, sampling of the signal represented by the function f(t) means an 
approximation to the function using a set of its values fk = f(tk) at the given set 
of the instants tf.. On the other hand, there is no reasonable way to deal with the 
integrals = / S(t — tk)f(t)dt using real sampler. We introduce a different method 
of discretization, based on a function exponentially localized both in the time and the 
frequency, instead of using Dirac delta function. The method is a realization of the 
representation that uses Gabor's idea of 'elementary' signals, which helps to overcome 
the difficulties with classical Gabor decomposition, connected with the incompatibility 
of orthogonality and localization, simultaneously, according to the Balian-Low theorem. 
Close to our approach is the one using Wilson bases [T3] . 

An introduction to the method of asymptotically exact expansions in TF plane with 
strong exponential localization property is provided below, followed by its applications 
to signal analysis, in a form of several basic examples, and using graphical representation 
constructed. 

2. Mathematical Preliminaries: the Unique Waveform and Asymptotically 
Exact Expansions 

An insight into the noncommutative geometry of the TF plane, as viewed from the 
position of signal processing, may be found by taking into account the fact that 
translations and modulations form a non-commutative continuous group: the Weyl- 
Heisenberg group. Its maximal commutative subgroup 'forming' the von Neumann 
lattice, called also Gabor frame, is subject to the restrictions imposed by the Balian- 
Low theorem. Consideration of the anticommutative subgroup of the Weyl-Heisenberg 
group, with double- density lattice, shows a way to deal with the singularities present. 
The subgroup under consideration is generated by translations, T a : f(t) — > f(t + a), 
and modulations, T b : f(t) — > exp(ibt) f (t) , with a = b = 1/7?, and anticommutation 
relation 



holds. The details of this approach are given in [TTJ. Because of the incompatibility 
of localization and orthogonality, one has either to use the description in terms 



T a T b + T b T c 
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of nonorthogonal subspaces, or to tolerate very weak localization (as for the 'sine' 
function), and infinite uncertainty products AuAt. 

The overall complication thus reduces to the consideration of four mutually 
nonorthogonal subspaces of the Hilbert space, represented by the respective four 
sublattices in the TF plane. Within each of the four subspaces there exists an 
orthonormal basis with strong exponential localization property, represented by the 
shifts of a single waveform, 



oc 



2 £ (-l) n a n e-™- 2 ™ 2 cosh (2/1*^) 

a(t) oc El . ( 2 ) 

eV^ 3 (t0F) 

The waveform is basic to our construction. The factor e~ nn in the sum above is 
responsible for strong exponential localization. The function (J2J) coincides with its own 
Fourier transform. (The same property has Hermite functions under proper scale, and 
that symmetry is exploited in our approach, see the next section.) The four sublattices 
correspond to the subspaces C H, fi = 1, 2, 3, 4, which i) are invariant with respect 
to T 2 and T b 2 , ii) M 3 has an orthonormal basis generated from the base waveform by the 
all combinations of the shifts T 2m and T b 2n with integer m and n (the third sublattice), 
and iii) M 4 = T a M 3 , M 2 = T b M 3 , and M 1 = T a M 2 = T b M 4 . 

The function #3 (a;) that enters the definition above corresponds to the so called 
lemniscate case of the theory of elliptic functions, and by #3 (a;) we mean Jacobi theta 
function ^ 3 (x, e _7r ); explicitly, 

00 

■& 3 (x) = l + 2^e -7m2 cos(2na;). (3) 

n=l 

The coefficients a n are defined as coefficients of the Fourier series expansion for the 
inverse square root of Jacobi theta function $4(0;) = $3(0; + tt/2): 
1 00 

oc 1 + 2 Y a n e~™ cos(2nz) . (4) 

The first six values of a n are: {0.501 052, 0.375 586, 0.312 961, 0.273 833, 0.246 446, 
0.225 907}. It may be said that the numbers are structure constants of the TF plane 
with respect to the symmetry group of that plane, the Weyl-Heisenberg group, that are 
generated by its anticommutative subgroup. 

For practical calculations, the function (J2J) can be rewritten in a form 

e~T + j: {-l) n a n e-™{e-- — ^ + e~- — ^) 
a(t) oc ^ . . (5) 



Comparison of the function (JSJ) with the Gaussian waveform is shown in FigJTJ Also 
shown are the three- and the five-term approximations to the function fl5J), in order to 
demonstrate the super convergence of the sum in the nominator. 
The phases of the basis vectors are fixed as follows: 

a m , n {t) = e inb{t+ma/2) a(t + ma), m, n = 0, ±1, ... . (6) 



Asymptotically Exact Localized Expansions for Signals in Time- Frequency Domain 4 




Figure 1. Power attenuation-both with the time and the frequency, for the signal ([5]), 
as compared to that for the Gaussian signal (the parabola) . Also shown are the three- 
and the five-term approximations to the signal (dashed lines). 



The most important feature clearly seen in the FigJT] is the strong exponential 
decay. The fact of the linear relation between the time and the frequency scales with 
attenuation is basic to our asymptotically exact reconstruction method, which follows 
from the resolution of identity for the Hilbert space of signals with finite energy. The 
respective scalar products are calculated as integrals 

f m ,n = J f(t)a m>n (t)dt, m,n = 0, ±1, ... . (7) 

The same fundamental fact of strong exponential localization allows to use finite 
sums as asymptotically exact approximations, 

/(*) ~ jE/^vWi (8) 

within the rectangle \m\ < M, \n\ < N. The finite sets {f m ,n} would represent the 
asymptotic expansion, understood in the sense of Hilbert space metrics. 

Because of nonorthogonality between the four sublattices, the energy of the signal 
is counted effectively twice so that 

\Y.\fmA 2 = J \f(t)\ 2 dt } • (9) 

With the equations ((?]) and flB]), it doesn't matter which particular representation of 
the signal under consideration is used, either the time or the frequency one. The numbers 
M and N are chosen such that the energy of the signal is well approximated according 
to the ()9]) with respect to the noise level. Further details of the basis construction, and 
its properties, can be found in the two papers [H] and [12] . 

Hilbert space vectors corresponding to the signals, in the basis under discussion are 
conveniently represented by complex matrices indexed by discrete time and frequency m 
and n, so that the corresponding "matrices" themselves play the role of the coordinates 
of the signal in the functional space. The "matrices" based on expansion formulae ([7]) 
and (IE]) can be represented in a graphical form, by using polar decomposition of complex 
numbers-the elements of the matrix, with semicircles positioned in the proper cells in 
TF plane, as described below. 
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Figure 2. Graphical representation for the function a(t) according to ([7]) and (|8]l. 
Time and frequency axes are horizontal and vertical, respectively. 

3. Graphical Representation: Sampling and Reconstruction of Commonly 
Used Signals 

Graphical representation for the function ([5]) itself that is based on asymptotic expansion 
of the type ([8]) is shown in FigEJ This example is basic: it makes explicit the 
noncommutative structure of the plane that is reflected in the nonorthogonality between 
the four subspaces of the Hilbert space of signals, which are represented by the 
respective four sublattices marked by alternating levels of grayscale. Recall that the 
symmetry of the four sublattices is described by the anticommutative pair of T a and 
Tb (DQ). The expansion is characterized by the complex numbers-the amplitudes (|7|), 
that are attributed to each site of this double overcomplete union of sublattices. 
The interpretation of the graphs introduced is straightforward: the sites of the four 
sublattices represent bases, orthonormal within each of the four subspaces. The sites 
with coordinates {ma,nb} are attributed to the functions a m ^ n (t) ([6]). The radii of 
the circles are equal to the absolute values \fm, n \, and rotations represent their phases, 
Arg(f m ^ n ); black semicircle down means zero phase. For instance, in Figj2]the amplitude 
at the origin, /o,o ; is real and positive. The nonzero amplitudes on the adjacent lattices 
reflect the nonorthogonality of the four subspaces of the Hilbert space of signals, which 
is inevitable in the view of the noncommutativity of the TF plane. Considered here 
are expansions for normalized signals (signals of unit energy); the portion of the energy 
of the signal at the given site is proportional to the respective circle area. Symmetry 
with respect to Fourier transformation is reflected in the symmetry with respect to the 
time and the frequency axes. Graphical representation described here demonstrates the 
functionality of our method based on expansion formulae (I7|) and ([8]). 

A convenient tool for the representation of common waveforms are the expansions 
in Hermite functions; the latter form orthonormal basis in the space of signals with 
finite energy. On the other hand, it is advantageous to use their asymptotic expansions 
of the type flB]) for the representations of the waveforms translated-modulated. The 
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Figure 3. Gaussian waveform represented by its asymptotic expansion in the TF 
plane. 
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Figure 4. The difference between the base waveform ([5]) and the Gaussian waveform, 
both normalized. The radii of the circles are shown with magnification 30, to expose 
minor differences that cannot be seen clearly in Fig. 2 and Fig. 3. 



interrelation between the two bases, the one, represented by the four sublattices, and the 
other, represented by the set of Hermite functions, once established, offers a convenient 
tool for signal representation and analysis: the advantages of the latter with respect 
to the signals 'centered' about the origin now can be combined with the translation 
invariance of the former, both with the time and the frequency. 

Signals most localized in the TF plane are represented by the Gaussian waveform: 
its uniqueness is guaranteed by the minimization of the Heisenberg uncertainty relation. 
Less localized are the so-called Gaussian monocycles, and the pulses of higher order, 
represented in terms of quantum mechanics as higher excited states of a harmonic 
oscillator, and their hierarchy constitutes a good classification scheme for the degree 
of localization. The character of localization of the Gaussian signal is seen clearly in 
the Figj3j The main part of the energy of the signal is concentrated on the central sites 
of the lattice, and the other amplitudes nearly vanish. The rotational symmetry of the 
expansion is an effect of the symmetry with respect to the Fourier transformation. The 
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Figure 5. Gaussian monocycle waveform in the TF plane. 

position and labeling of the sublattices coincide with those in the Figf2j The minor 
differences between a(t) and the Gaussian waveform that reflect the orthonormality of 
the former on its own sublattice, and the absence of such for the latter, are plotted 
in FigJH Notice that the difference maxima are on the sites (±2a, 0), and (0, ±26), 
which belong to the 'central' sublattice, and that the central difference is relatively very 
small. The next waveform considered is the so called Gaussian monocycle, FigfSJ It 
is less localized than the Gaussian one, as expected. The symmetry of the figure is 
due to the symmetry of the waveform with respect to the Fourier transformation, as 
for the Gaussian waveform. In contrast with the Gaussian one, the central amplitude 
is exactly zero, and the amplitudes on the frequency axis are purely imaginary; the 
phases of the amplitudes on the other sites vary in a symmetric manner. The two 
waveforms considered above are actually the zeroth and the first in the set of Hermite 
functions, which constitute a basis in the Hilbert space of signals, as it was mentioned 
above. Since the interrelation of the two complete sets is of primary interest, the other 
Hermite functions will be considered too. The hierarchy of the eigenvectors of quantum 
harmonic oscillator Hamiltonian can represent the degree of localization, parameterized 
by the order of the polynomial, /. Normalized Hermite functions are 

u,(t) = - T ^=e~ t2 / 2 H l (t), 1 = 0,1,2,... (10) 
V^!2V^ 

and Hi(t) are Hermite polynomials 

m) = (-l) l e t2 ^e- t2 . (11) 

The main features of the graphical representation for the two waveforms considered 
above, their symmetry and localization, are persistent with the whole hierarchy. Hermite 
functions of higher orders are shown in Figj6j By establishing the interrelation between 
the two bases, the one, localized about the origin, and the other, uniformly spaced in 
the plane, we show the way to localized asymptotic approximations applicable to any 
signal. 
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Figure 6. 


Hcrmitc functions of higher orders 


n TF plane. 









With the property of the strong exponential localization of the waveform (jSJ), the 
expansion method presented is equally applicable to arbitrary displacements in the 
TF plane. In many signal processing applications the recognition of a standard signal 
displaced in the TF plane, i.e. translated and modulated, is sought. Calling it a(t), the 
problem reduces to calculations with single function of the displacement, C(£, r/), that 
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Figure 7. Standard signal at a 'random' position in TF plane so close to the origin 
that the uncertainty product is much less than the uncertainty limit imposed by the 
Heisenberg relation. The signal expansion is shifted intentionally according to the 
values of the At and Aw. 

is represented by the integral 

C m , n (£, V) = J e^ t+ ^a(t + O^Jt)dt. (12) 

It is clear that the restrictions of the uncertainty relations would not affect the 
recognition, within our method using expansions of the type (JSJ). On the other hand, the 
set of scalar products represented by the integral (fT2"|) reveals noncommutative geometry 
of the TF plane in a manner different from that of the FigfSJ Fig J7] presents an example: 
the standard signal at 'random' position in the TF plane very close to the origin, which 
gives another insight into the structure of the space of signals in quantum uncertainty 
limit. 

4. Discussion and Conclusions 

The sample of the signal represented by the function f(t), when understood directly, 
is the value of the function at the given instant to- Clearly, such understanding is 
not an 'operational' one. By assuming the form of integrals (J7J) for the samples, the 
concept of sampling in TF plane looks much more transparent; and it is operational. 
In this course, the conceptual aspect may be formulated simply as a rule to deal 
with operational quantities. Furthermore, such an approach is a realization of the 
representation that uses Gabor's 'elementary' signals, which helps to overcome the 
difficulties with classical Gabor decomposition, connected with the incompatibility of 
orthogonality and localization, simultaneously, according to the Balian-Low theorem. 
The inevitable complications caused by the partial nonorthogonality then reflect the 
noncommutativity features of the TF plane. It should be stressed that the scalar 
products represented by (0) have invariant meaning. 
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In conclusion, (i) an asymptotically exact mathematical method for localized 
expansions of signals in TF domain is presented; (ii) based on the method, a graphical 
representation for signals in the TF plane is introduced; (iii) noncommutative geometry 
of the TF plane is visualized; (iv) the interrelation with the Hermitian basis is 
established, in terms of the graphical representation constructed. 
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